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Abstract. There has recently been considerable interest in quantum walks in connection with 
quantum computing. The walk can be considered as a quantum version of the so-called correlated 
random walk. We clarify a strong structural similarity between both walks and study limit theorems 
\q [ and absorption problems for correlated random walks by our PQRS method, which was used in 
our analysis of quantum walks. 
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The correlated random walk in one dimension can be denned as follows. We suppose that a 
particle performs a random walk on a state space 5*. In this paper, we consider three cases 
as 5*, that is, S = Z, S = Z + , and S = {0, 1, ... , iV}, where Z is the set of integers and 
Ch ' Z + = {0, 1, . . .}. The evolution of the correlated random walk is given by 

q-( P{ the particle moves one unit to the left) 

p, if the previous step was to the left, 
^ ; [ 1 — q, if the previous step was to the right, 

H 

and 

P( the particle moves one unit to the right) 

1 — p, if the previous step was to the left, 

q, if the previous step was to the right. 

When p = q, the probability that the particle moves one unit in the same direction as the 
last step is p, and the probability that the particle moves one unit in the opposite direction 
as the last step is 1 — p. In particular, if p — q — 1/2, then the walk is equivalent to the 
well-known symmetric (non-correlated) random walk, i.e., the particle moves at each step 



Abbr. title: Limit theorems and absorption problems for correlated random walks. 
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either one unit to the left with probability 1/2 or one unit to the right with probability 1/2. 
The directions of different steps are independent of each other. 

The correlated random walk has been studied by various authors. Recent work relating 
to the present paper includes [21 El [Zl El El HOl EOl E21 E3l [26] - More detailed information on 
their results will be stated in the next section. 

On the other hand, a discrete-time quantum walk was introduced and intensively studied 
by [U E] as a quantum variant of the classical random walk and since then considerable 
work has been done on quantum walks by a number of groups in connection with quantum 
computing (see [H [131 ES], for a review). Our group has also been studying quantum walks 
by using mainly the PQRS method based on combinatorics which is different from Fourier 
analysis, (see Konno [HI H21 HB] , Konno, Namiki and Soshi [17J, Konno et al. [18J). Meyer 
and Blumer [21] pointed out that the quantum walk can be considered as a quantum version 
of the correlated random walk. More recently, Romanelli et al. [24J analyzed the quantum 
walk by separating the quantum evolution into Markovian and interference terms. They 
reported that if the interference terms are neglected then the resulting equation becomes 
the Telegraphist's equation which can be considered as a space-time continuum limit of the 
correlated random walk. It is a well established result for a one-dimensional quantum walk 
that if (i) it is expressed in terms of pairs of real probabilities and (ii) quantum correlations 
are neglected, then, it reduces to a two-dimensional Markovian process which is equivalent 
to a one- dimensional correlated random walk (or, equivalently, to a Telegraphist's equation 
in the continuum limit). The same fact holds for a one-dimensional correlated (or persistent) 
random walk: it is equivalent to a Markovian process in two dimensions. As the following 
definitions show, there is a strong structural similarity between them. One of our motivations 
for this study is to clarify the gap between them by using our PQRS method. We present 
limit theorems and absorption probabilities for correlated random walks as in the case of 
quantum walks (see [16]). 

This paper is organized as follows. Section 2 treats the definitions of correlated random 
walk and the quantum walk, and explains the PQRS method. Moreover we summarize some 
recent work on correlated random walks. In Section 3, we give two types of limit theorems. 
Section 4 is devoted to one- and two-boundary absorption problems. 



2 Definition and PQRS Method 

First we present the definition of the correlated random walk, which clarifies a connection 
with the quantum walk. The time evolution of the walk is given by the following transposed 
transition matrix: 



A 



a b 
c d 



where 



a = p, 6 = 1 — q, c — 1 — p, d = q. 

In order to avoid trivialities, we assume that < a, d < 1 throughout this paper. In 
the case of quantum walks, A is a unitary matrix with a, b, c, d e C, where C is the set of 
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complex numbers. By the unitarity, elements of A satisfy the following relations: \a\ 2 + \c\ 2 



\b\ 2 + \d\ 2 = 1, ac + bd 



0, c 



-Ab, d = Aa, where z is the complex conjugate of z G C 



and A = det A = ad — be. 

For both correlated random walks and quantum walks, we can think of the evolution as 
driven by the outcome of tosses of two coins, with each coin determining the evolution at the 
time step for which that coin is active. The directions R and L can be respectively thought 
of as corresponding to the heads and tails of the coin, or equivalently as an internal chirality 
state of the particle. The value of the coin controls the direction in which the particle moves. 
In fact, when the coin shows L, the particle moves one unit to the left, when it shows R, it 
moves one unit to the right. Let 
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R 




1 



so we have 



UL = aL + cR, 
UR = bL + dR. 

From now on, in order to emphasize the connection with the quantum walk, we sometimes 
use notation taken from this field. Let \I/fc(n)(£ R 2 ) be the probability that the correlated 
random walk is in location k at time n, 



**(n) 



with the chirality being left (upper component) or right (lower component). Then the 
dynamics for \l/fc(n) is given by the following transformation: 



^ k (n + 1) = P V k +i{n) + Q * fc _i(n) 



(2.1) 



where 



P 



a b 




Q 





c d 



It is noted that A = P + Q. Romanelli et al. [21] gave a similar argument, so we clarify 
the relation between our setting and theirs. In our notation, the model they considered 
corresponds to the following correlated random walk: 



A 



cos 2 9 sin 2 6 
sin 2 9 cos 2 9 



that is, a = d = cos 2 9 and b = c = sin 2 9 for some value of 9. Therefore, by using Eq. (12. ip . 
we have 



^(n + 1) = cos 2 9 *£ +1 (n) + sin 2 9 
*(n + 1) = sin 2 9 + c os 2 9 *J_ x i 



n 



n 
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The above equations can be also obtained when the interference terms are neglected in 
Eq. (5) in their paper. We remark that ^ 3 k {n) (our notation) is equivalent to Pkj(n) (their 
notation) (j = L,R). 

One of the essential parts of the present paper is the study on the dependence of some 
important properties and quantities (e.g., limit distribution, symmetry of distribution, ab- 
sorption probability) on the initial distribution. Therefore we introduce the collection of 
initial distributions as follows: 



: a + p = 1, a,0 > 



Let X% be the correlated random walk at time n starting from initial distribution tp e $ 
with Xq = 0. Let S(/,m) be the sum over possible paths for which the particle arrives at 
location k at time n starting from the origin with I + m = n and —l + m = k. For example, 

5(3, 1) = QPPP + PQPP + PPQP + PPPQ. 

Then 





a 


{'- 





P{X* = k) = [l,l] S(J,m) 



a 



where 



x. 



l [u,v\) = xu + yv, 1 = *[1, 1], ip — *[ct, j3] G R 2 is an initial distribution with 
a + (3 = l,a, (3 > 0, and t means the transposed operator. On the other hand, in the case 
of quantum walks, 



P(XZ = k) = \\E(!,m)<p\\ 



where ip = t [a,f3] e C 2 is the initial qubit state with \a\ 2 + \f3\ 2 = 1 and 



X, 



\y\ 



In our treatment of correlated random walks, as well as the matrices P and Q, it is 
convenient to introduce 

R 

When ad — be (= a + d — 1) ^ 0, P, Q, R, and S form a basis of M 2 (R) which is the vector 
space of real 2x2 matrices. Now, we let 



c d 




"00" 





, s = 


a b 



X 



x y 
z w 



If ad — bey^ 0, then 

X = c p P + c q Q + c r R + c s S, 
where the coefficients c p , c q , c r , c s are determined by 



Cp 




d 


-c 




X 




1 


-b 


a 




y 


C s 


ad — be 





d -c 




z 


. C 1 . 







-b a 




w 
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Therefore we can express any 2x2 matrix X conveniently in the form, 

X = 1 {{dx - cy)P + (-bz + aw)Q + (-6a; + ay)R + (dz - cw)S} . (2.2) 

CLCL C/C- 

From now on, we assume that ad — be ^ 0. We call the analysis based on P,Q,R, S the 
PQRS method in this paper. 

The n x n unit and zero matrices are written as I n and O n respectively. If X = I 2 , then 
Eq. (T22D gives 

! (dP + aQ-bR-cS). (2.3) 



ad — be 

The following table of products of P,Q,R, S is very useful in computing some quantities 





P 


Q 


R 


S 


p 


aP 


bR 


aR 


bP 


Q 


cS 


dQ 


cQ 


dS 


R 


cP 


dR 


cR 


dP 


S 


aS 


bQ 


aQ 


bS 



where PQ = bR, for example. 

In order to consider absorption problems in Section 4, we describe the absorption prob- 
ability of a correlated random walk starting from location k on {0, 1, . . . , iV} (N < oo) with 
absorbing boundaries and iV as follows. 

First we consider N < oo case. Let (n) be the sum over possible paths for which the 
particle first hits at time n starting from k before it arrives at N. For example, 

2! 3) (5) = P 2 QPQ = ab 2 cR. 

The probability that the particle first hits at time n starting from k with initial distribution 
<p(£ $) before it arrives at iV is 



P ( k N \n^) = \H N \nM 



l- 



So the probability that the particle first hits starting from k with initial distribution ip(E $) 
before it arrives at N is 



n=0 

Next we consider N = oo case. Let E$°\n) be the sum over possible paths for which 
the particle hits at time n starting from k. For example, 

Hj 00 ^) = P 2 QPQ + P 3 Q 2 = (ab 2 c + a 2 bd)R. 

In a similar way, we can define Pjf°\n; tp) and Pjf°\tp). 

From now on we review some recent works on correlated random walks. Renshaw and 
Henderson (23] considered the a = d case and obtained an expression for P(X% = k) and its 
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limiting distribution as n — > oo with a — > 1 for the state space S = Z. Zhang [26J investigated 
the absorption probability and the expected duration for various types of boundaries such as 
an absorbing boundary or an elastic boundary. Lai and Bhat [2D] calculated P(X% = k) for 
S — Z, and limiting distribution for S = {0, 1, . . .} and S — {0, 1, . . . N}. Bohm [8 J consid- 
ered the a = b case and studied the distribution of the absorption time and the asymptotic 
distribution of the maximum for this walk on S = {0, 1, . . .}. Bohm [7j computed P(X% = k) 
for S — Z, {0, 1, . . .}, and {0, 1, . . . N} based on Krattenthaler's result (Krattenthaler [19J) 
for counting lattice paths with turns. Moreover, when a = d, Bohm gave the asymptotic 
behaviour of P[X^ = k) as n — > oo. In addition, Allaart and Monticino [3] analyzed optimal 
stopping rules for a general class of correlated random walks. More refined results were 
obtained by Allaart |2J. Chen and Renshaw investigated correlated random walks in higher 
dimensions jH] and considered a more general walk [TUJ. Mukherjea and Steele [22J studied 
a correlated gambler's ruin problem. 

3 Characteristic Function and Limit Theorem 

This section treats two different types of limit theorems for correlated random walks X%. To 
do so, first we need to know P(X% = k) for n + k = even. For fixed / and m with I + m = n 
and — I + m = k, we have 

S(Z,m)= Yl P h Q mi P h Q m2 ■■■P l "Q mn , 

/j,m;>0:rrtiH \-m n =m,l\-\ \-ln=l 

since 

?(ir = fc)=*lS(!,m) ip. 
Noting that P, Q, R, and S are a basis of M 2 (R), 2(7, m) has the following form: 

S(7, m) = p(l, m)P + m)Q + r(Z, m)R + s(l, m)S. 

Our next aim is to obtain explicit forms for p(l, m), q(l, m),r(l, m), and s(l, m). In the case 
of n = I + m = 4, we know 

5(4, 0) = a 3 P, 3(3, 1) = 2abcP + a 2 bR + a 2 cS, 
3(2, 2) = bcdP + abcQ + b(ad + bc)R + c(ad + bc)S, 
3(1, 3) = 2bcdQ + bd 2 R + cd 2 S, 3(2, 2) = d 3 Q. 

So, for example, 

p{3,l) = 2abc, g(3,l) = 0, r(3, 1) = a 2 b, s(3, 1) = a 2 c. 

Now the following key lemma can be obtained by using our combinatorial method. 

LEMMA 1 We consider correlated random walks with abed ^ 0. Suppose that l,m > 
with I + m = n, then we have 
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(i) for I A m(= min{Z, m}) > 1, 



(ii) for l(= n) > 1, m = 0, 
fmj /or / = 0, m(= n) > I, 

The proofs of parts (ii) and (iii) are trivial. The proof of part (i) is based on enumerating 
the paths of drift I + m = n. To do so, we consider the following 4 cases: 





U'l 


w 2 


W3 




W2 1 






p(l, m) : 


pp- 


~PQQ^ 


~QPP^ 


~P- 




~QPP^~- 






Wl 


W2 








W2~/+l 




q(l,m) : 


QQ^ 


~QPP^ 


~PQQ^ 


~Q- 


P~P^~ 


~PQQ^~- 






Wl 




W3 










r(l, m) : 


p~p^~ 


~P / QQ^ 


~^PP^ 


~P- 












LL'2 












s(l, m) : 




~$pp^ 


~P'QQ^ 




P~P^~ 


~P, 





where Wi, ty 2 , • • • , W27+1 — 1 an d 7 > 1- We remark that p(l,m),q(l,m),r(l,m), s(l,m) in 
S(Z,m) correspond to the types of paths described above in terms of P and Q, respectively. 
A similar argument can be seen in the proof of Theorem 2.1 of Bohm [7j. On the other hand, 
for the quantum walk case, this proof appears in Konno [15]. One of the main difference 
between the correlated random walk and quantum walk is the sign of the term (be/ ad). That 
is, 

be _ ( (1 — a)(l — d)/ad > 0, correlated random walk, 
ad \ —(1 — |<i| 2 )/|a| 2 < 0, quantum walk. 

The negativity in the quantum walk case corresponds to the interference effect. For the 
convenience of readers, we give a outline of the proof based on [15]. 

Proof. Here we consider only p(l,m) case. Since the proofs of the other cases are similar, 
we will omit them. First we assume I > 2 and m > 1. In order to compute p(l,m), it is 
sufficient to consider only the following case: 

Wl W 2 UI3 U>2~f W27+I 

C(P, w) ( ^ +1) (l, m) = P~P^QQ^~QP~P^~P- ■ ■ QQ^~QPP^~P, 

where w = (wi, W2, ■ ■ ■ , W2 7 +i) Z+ 7+1 with w%, w 2 , ■ ■ ■ , W27+1 > 1 and 7 > 1. For example, 
the string PQP has w\ = w 2 = W3 = 1 and 7 = 1. We should remark that / is the number 
of Ps and m is the number of Qs, so we have 

I = Wi + W 3 H h U>2 7 +1, 

m = w 2 + W4 H h w 2l . 



07 0(7 c 



S(Z,0) = a 1 ' 1 P. 
S(0,m) = d m - l Q. 
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Moreover 27 + 1 is the number of clusters of P's and Q's. Next we consider the range of 7. 
The minimum is 7 = 1, that is, 3 clusters. This case is 

P...PQ...QP...P. 

The maximum is 7 = (/ — 1) A m. This case includes the patterns: 

PQPQPQ ■ ■ ■ PQPQP ■ ■ ■ P (I - 1 > m), PQPQPQ ■ ■ ■ PQPQ ■ ■ -QP (I - 1 < m), 

for examples. Here we introduce a set of sequences with 27 + 1 components: for fixed 
7 G [1,(1-1) Am], 

W(P, 27 + I) = {w = (w 1: w 2 , • ■ ■ , w 2 y+i) G Z 27+1 : u>i + u>3 H h ^27+1 = I, 

w 2 + w 4 H h W27 = W, lUi, w 2: . . w 2y , w 2l+ i > 1}. 

From the following relations: 

P 2 = aP, Q 2 = dQ, PQ = bR, R 2 = cR, RP = cP, 

we have 

C{P,w)^ +1) (l,m) = a Wl - 1 Pd W2 - 1 Qa w *- 1 P---d w ^- 1 Qa W2 ^ 1 - 1 P 
= a l -^ +1) d m -^{PQyP 
= a l -^ +1) d m ~WR?P 
= a l -^ +1) d m -Wc<- l RP 
= a l -^ +1) d m -Wc<P, 

where w G W(P, 27 + 1). For I > 2, m > 1, that is, 7 > 1, we obtain 

C(P,w)^ +1 \l,m) = a'- (7+1) 6 7 c 7 d m - 7 P. 

Note that the right hand side of the above equation does not depend on w G W(P, 27 + 1). 
Finally we compute the number of w — (wi,w 2 , . . . , w 2l , u>2 7 +i) satisfying w G W(P, 27 + 1) 
by a standard combinatorial argument as follows: 



\W(P,2 1+ 1) 



From the above observation, we obtain 



Z - l\/m - 1 
7 / \7 - 1 



(i-l)Am 

p(l,m)P = J2 E C(P,w)^ +1 \l,m) 

7=1 w&W{P,2-y+l) 
(i-l)Am 

= £ l^2 7 + l)|C(P, W )i 27+1 )(/,m) 
7=1 

(i-l)Am 



E ( z 7 x ) (™_ 1) «^ (7+1) & 7 c 7 rf m - 7 p. 



7=1 
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So we conclude that 

p(l,m) 

When / > 1 and m = 0, it is easy to see that 



p(l,0)P = P l = a 1 ' 1 P. 
Furthermore, when I = 1, m > 1 and / = 0, m > 0, it is clear that 

p(l, m) = 0. 

Similarly we have the following results for q(l,m),r(l,m) and s(l,m): if I > 1 and m > 2, 
then 

g(i,m)= £ (^ J( m ~ ) ai_Wrfm " (7+1) ' 

if / = and m > 1, then g(0, m) = d m_1 . When I > l,m = 1 and / > 0, m = 0, we know 
q(l,m) =0. If Z,m > 1, then 

/_1 V rn_1 V^ 



if I A m = 0, then r(l, m) = 0. If I, m > 1, then 

.(/,„,) =E(': 1 1 )(;: 1 1 ) o '" 767 " V7dro " 7 ' 

if / A m = 0, then s(7, m) = 0. Therefore, when I A to > 1, we have 



lAm 



E(l,m) = a l d rn J2 



bc\ ( I — 1 \ /m — 1 



-'Kadi \7 — l/\7 — 1 

7=1 v ' v / \ ' 

So the proof of Lemma [T] is complete. 
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In order to obtain an expression for the distribution of X^, we will use the following 
facts: for k = 1, 2, ... , [n/2], 

V /&cV _1 1 fn-k-1 



/ 6c\ 7 1 fk — \\ fn — k — 1 



^ ( ad 



7=1 



7 — 1/ V 7 — 1 



., J = 2 i 71 i(-(A; - 1), -{(n - A;) - 1}; 2; be/ 'ad), 
2 Fi{-(k - 1), -{(n ~k) - 1}; 1; 6c/ad), 



where [x] is the integer part of x and 2 i*i (a, o; c; z) is the hypergeometric series. In general, 
as for orthogonal polynomials, see Andrews, Askey and Roy [6], for example. Then the 
distribution of X% can be derived from Lemma [1] as follows. 
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LEMMA 2 For n > 2 and k = 1, 2, ... , [n/2] 

P(X* = n-2k) 
= a k ~ 2 d n - k - 2 



bc{(dk + c(n - k))aoc + (Wfe + a(n - A;))d/3jF 2 (ri ' fc) 

+ (aca + bd(3)AFi n ' k) 



P{X* = -{n-2k)) 



k-2m-k-2 



6c|(c/c + d(n — fc))ao! + (a/c + 6(n — k))d(3^F : 



(n,k) 



for n > 1, 



+ (aca + bd(3)AF{ 



P(X% = n) = d n -\ca + d/3), P(X% = -n) = a n ~ l (aa + b/3) 



(n.k) 



where A = ad — be, and F i 



(n,k) 



iFi(-(k - 1), -{{n - k) - 1}; i; be/ ad) (i = 1,2). 



The above result for ip = *[0, 1] and </? = *[1, 0] corresponds to Eqs. (19) and (20) in Corol- 
lary 2.6 of [7], respectively. By using Lemma HJ we have an expression for the characteristic 
function of Xf. This result will be used in order to obtain limit theorems of Xf. 



THEOREM 3 When abed ^ 0, we have 



E( e * x " 



a n ~\aa + b/3) + d n -\ca + dp) 
+i 



-a n -\aa + b(3) + d n -\ca + 
+ £ a fc - 2 d 



cos(n£) 
sin(n£) 



jn— fc— 2 



fe=l 



i(n,k) 



+i 



( n 


n 








■•) 


\2 " 


2_ 



6cn| (c + d)aa + (a + 5)<i/3 j.F 2 

+2(aca + &d/3)AF 1 (n,fc) l cos((n - 2Jfe)f) 
6c(n - 2A;){(c - d)aa + (a - 6)d/3}F 2 ( "- fc) sin((n - 2k)£) 

n/2-2 



ben 



{(c + d)aa + (a + 6)d/3} F 2 (n,fc) + (aca + bd(3)AFl 



i(n,k) 



where l(x,0) = l(resp. = 0) if x = (Vesp. x ^0). 

From this theorem, we have the mth moment of X% in the standard fashion. The following 
result can be used in order to study symmetry of distribution of X%. 
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COROLLARY 4 We assume that abed ^ 0. 
(i) When m is odd, we have 



ZD 



a n ^(aa + b/3) + d n ~\ca + dp) 



n 



2 1 



+ a k - 2 d n ~ k ~ 2 bc{n-2k) m+l ^(c-d)aa 



[a — 



k=l 



b)d/3^F^ h) . 



(ii) When m is even, we have 



tp\m\ 



.n-1 



(aa + b/3) + d n ~^(ca + d(3) 



n 



+ J2 ak ' 2d 



k-2 m-k-2, 



n - 2k) 1 



fc=l 



ocn{(c + d)aa + (a + b)d(3^F^ n ' k) + 2(aca + bd(5)/\F[ 



(n,k) 



By using Lemma [2] and CorollayH] (1) with m = 1, we obtain the following necessary and 
sufficient condition for symmetry of distribution of X% with a = d. 

PROPOSITION 5 (i)Ifa = de (0, 1) and a ^ 1/2, then 

$, = $ = {^ = t [l/2,l/2]}, 

(ii) if a = d = 1/2, then 



$ = $, 



where 



$ 



{(p E $ : P{X% = k) = P(X% = -k) for anyn E Z + andk E Z}, 
{ip E $ : £7(X£) = /or anyn E Z + } . 



The dependence of the symmetry of the distribution on the initial condition is remarkably 
different between classical and quantum walks. For results on the symmetry of the quantum 
walk, see [TJJ [El EH E] • I n the last part of this section, we will treat two types of limit 
theorem for X% with different scalings. The first limit theorem for the case of X% / ^Jn shows 
that the correlated random walk considered here behaves like as a classical (non-correlated) 
random walk. The following result can be derived from a consequence of Theorem 3.1 in 
Bohm [7] with m, r, t — > 0, a — > a and z — > (1 — a) /a. 

THEOREM 6 Consider the correlated random walk with 

a 1 — a 
1 — a a 
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If n — > oo, then 



where the distribution ofW^ is N(0,a/(1 — a)), with N(m,cr 2 ) the normal distribution with 
mean m and variance a 2 and Y n =>- Y means that Y n converges weakly to a limit Y . 

It is noted that the above limit distribution is independent of the initial distribution 
ip £ Theorem corresponds to Eq. (16) in [23], since variance a/(l — a) becomes 
cos 2 9/(1 — cos 2 9) = cot 2 9 = D(9) in their setting. 

Next we consider how we renormalize X% when a — a n — > 1, that is, the above variance 
a n /(l — a n ) — > oo. The answer for this question is the following theorem. 

THEOREM 7 Consider the correlated random walk with 

Clrt. 1 CLn 



A, 



1 — CL r . 



a, 



where 



9 

1--, 

n 



(3.4) 



and a fixed 9 G (0, 1). Ifn—> oo, then 



_n ^ 

n 

where the probability measure of is the sum of its atomic part \i\ and the absolutely 
continuous part fi 2 as follows. The atomic part is 

where 5 X is the Dirac measure at location x. The absolutely continuous part \x 2 is given by 
the density function: 



i (9Vi - x 2 ) + 



i 



X/ 



= Ji(0VT^x2) 



for x G (—1, 1), where I v (z) is the modified Bessel function of the order v, that is, 

(z/2) v+2n 



n=0 



n\ T(u + n + l) 



and T(z) is the gamma function. 

Proof. We begin with an asymptotic result on the hypergeometric series: 

LEMMA 8 Let a n = 1 - (6/n) with < 9 < 1. If n oo with k/n = x G (0, 1/2), then 



2 F 1 (-(k - 1), -{(n -k)- 1}; 1; (1 - a n ) 2 /a 2 n ) -> / O (20v / T^ 2 ) 



2 F 1 (-(k - 1), -{{n - k) - 1}; 2; (1 - a n ) 2 /a 2 n ) 



9y/x(l — X) 



h^VT^x 2 ). 
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One of the above derivations can be seen on page 411 of [23] . By using Theorem [3] and Lemma 
[HI we obtain the following asymptotics of characteristic function E(e l ^ Xn ^ n ): if n — > oo with 
k/n = x e (0,1/2), then 

e~ e cos£ + ie~ e {— a + ($) sin£ 

•1/2 



+0e" 



Therefore we have 



2I (29y/x(l-x)) + 



y/x(l - x) 



h{26^x{l-x)) cos((l - 2x)£)dx. 



lim E(e*^) = e- e cos£ + ie- e (-a + [3)sm£ 



6e 



-i 



/o(^vT 



1 



x- 



-.I\{9y/l — x 2 ) cos(x^)dx 



{aS-i(x) +p6 1 (x)) 




where 



g(x)S a (x)dx = g(a). 



^ x dx, 



So the desired conclusion is obtained. 



Interestingly, Theorem [7] shows that the limit distribution is a mixture of quantum and 
classical parts. That is, the atomic part fix, having two peaks, corresponds to a quantum walk 
(as for the limit theorem for the quantum walk, see [HE)]) and the absolutely continuous 
part f/,2, having a bell shape like a normal distribution, corresponds to a classical random 
walk. 



4 Absorption Problem 

In this section we consider absorption problems for correlated random walk on state spaces 
{0, 1, . . .} or {0, 1, ... , N} by using the PQRS method as in the case of the quantum walk 
(see pi]). 

Before we move to the correlated random walk, first we describe the classical (non- 
correlated) random walk on a finite set {0, 1, . . . , N} with two absorption barriers at locations 
and iV (see [HI [T2], for examples). The particle moves at each step either one unit to the 
left with probability p or one unit to the right with probability q = 1 — p until it hits one 
of the absorption barriers. The directions of different steps are independent of each other. 
The random walk starting from k G {0, 1, . . . , iV} at time n is denoted by here. Let 

T m = mm{n >0:Y*=m} 
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be the time of the first visit tome {0, 1, . . . , N}. Using the subscript k to indicate Yq = k, 
we let 

Pf } = P fc (T < T N ) 

be the probability that the particle hits starting from k before it arrives at N. The 
absorption problem is also known as the Gambler's ruin problem. To obtain Pf; , we use 
the difference equation 

pr=ppr-i+iPiti (4.5) 

for 1 < k < N — 1 with boundary conditions: 

P (A ° = 1, J#° = 0- (4-6) 

Our basic strategy is to consider a similar equation in the correlated random walk case and 
to apply our PQRS method to it. 

From now on we focus on 1 < k < N — 1 case for the correlated random walk. So we 
consider n > 1 case. Noting that {P, Q,R, S} is a basis of M 2 (R), (n) can be written 

as 

Slf°(n) = p[ N \n)P + q[ N \n)Q + rf \n)R + s^\n)S. 

From the definition of (n), it is easily shown that there exist only two types of paths, 
that is, P — > ■ ■ ■ — > P and Q —>••••—>■ P, since we consider only a hitting time to before 
it arrives at N. Therefore we see that q k \n) = sjj. (n) = (n > 1). 
We assume N > 3. Noting that the definition of (n), we have 



a; ' n = a 



n(n-l)P + 8Ui(n-l)Q, 



for 1 < A; < iV — 1. The above equation corresponds to the difference equation, i.e., Eq. (14. 5 j) 
for the classical random walk. A similar approach can be seen in [2U]- Then we have 

(N) , \ (JV) / . (AT) / n 

rf°(n) = ^{(n - 1) + dr^n - 1). 

Next we consider boundary conditions related to Eq. (14. 6p as in the case of the classical 
random walk. When k = N, 

pf ) (0;^) = 4 i 3^(0)^ = 0, 

for any ip 6 <&. So we take Sj^ (0) = 2 . In this case, Eq. (12. 2p gives 

P i N N \0)=rP(0) = 0. 

If k = 0, then 

PfW) = 4 l 2^(0)^ = 1, 
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for any pG$. So we choose 5^(0) = From Eq. (I2.3p . we have 



iT(o) 



d 



ad — 6c 



a<i — 6c 



Let 



Pk ( n ) 



Therefore we can formulate the absorption problems as follows: for n > 1 and 1 < k < N — 1, 



a c 




-[(n-l) + 





b d 



and for 1 < k < N, 



4 N \o) 



1 



ad — be 



d 



vl N \0) 








(4.7) 



Moreover, 








for n > 1. So the definition of Pf; ((p) gives 



(4.i 



n=l 



where <£> = *[a, /3] € $ and 

Pf(n; = (aa + b(3)p[ N \n) + (ca + d/?)rf \n). 
To solve Pj: N \(f), we introduce generating functions of p k N \n) and r£ (n) as follows: 



(4.9) 



(AO/ 



71=1 



n=l 



From Eqs. (14. 8 p and ( 14. 9p . we obtain the following lemma: 

LEMMA 9 We assume that < a, d < 1 wift c = 1 — a, b — 1 — d and ad — be 7^ 0. T/ien 
we have 



Pf >(n; *>) = (aa + WT(1) + (ca + 



(4.10) 
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By Eq. (14. 7p . we have 

pf\z) = atf&tz) + cz7f_\(z), (4.11) 

rf\z) = bz$t\(z) + dzrf + \(z). (4.12) 

To solve these, we remark that both pj^\z) and r^\z) satisfy the same recurrence: 



(Az + iy»\(z) + ati N \z) = 0, 



d?Zl(z)-^Az + ^y^(z) + a^\z)=0, 

where A = det A = ad — be. The characteristic equations with respect to the above recur- 
rences have the same root: if < a, d < 1, then 



A2 2 + l^ v /A¥-2(a(l + 6c)z 2 + l 

A± " X±{Z) ~ 2dz ' 

From now on we consider the a = d 1/2) case with iV = oo. We remark that the 
definition of (n) gives p^in) = (n > 2) and Pi (1) = 1. So we have p^ 00 "*^) = z. 
Moreover noting \imk-^oaPk\ z ) < oo, < A + < 1 < A_ (0 < z < 1) and using Eq. (14.111) . 
we obtain the explicit form 

vt\z)=z\X-\ ^\ z) = h±Z^^-\ 

where 

1 + (2a - l)z 2 t V( 2a ~ 55 - 2 W + (1 - a) 2 }* 2 + 1 
± 2a^ 
Therefore we have p^ (1) = ^^(l) = 1 for any > 0. By using these and Lemma [9], we 
obtain Pj: ((p) — 1 for a = d 1/2) case. On the other hand, in the case of a = d = 1/2, 
the well known result for symmetric classical random walk and Proposition (ii) give the 
same conclusion. So we obtain 

PROPOSITION 10 Assume that a = d e (0, 1). For any k > and initial distribution 
<p = t [a, {3] G $, we /iai>e 

P fc (oo V) = 1. (4.13) 

Next we consider a finite iV case. First we note that A + A_ = 1 since a = d. Then pji. 
and satisfy 

pf = A 2 A+ _1 + £? z A fc -\ 7i N \z) = C z X k + - N+1 + D.A™. (4.14) 

By using Eq. (14. 14j) and boundary conditions: p\ (z) = z and r^li(z) = 0, we see that 
pf\z) and rl N \z) satisfy 

A N \z) = (| + £,) A^ 1 + (| - E,) A*" 1 , (4.15) 
rf>(z) = C z (\ k + - N+1 - A^ 1 ). (4.16) 
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On the other hand, from Eq. (14. 12p and r^\(z) = 0, we have 

rf\z) = z{(l - a)^ N \z) + a^f\z)}, ?fl 2 (z) = (1 - a)z^\{z). (4.17) 

Combining Eqs. (14.151) and (14.161) with Eq. (I4.17P gives 

C 2 = (l-a)^ 2 (A^- 3 -A^- 3 ) (4.18) 
x {-(A^ 2 - \ N - 2 ) 2 + a^(A^ 2 - A W - 2 )(A^ 3 - A^ 3 ) 

+(A + -A_) 2 }^, 



E 



2(A + -A_)(A^ 3 -A N - 3 ) (4.19) 



2 (A^ 2 - A^- 2 ) 
x {-(A^- 2 - A^ 2 ) 2 + az(X^ 2 - A^- 2 )(A^- 3 - A w ~ 3 ) 

+(A + -A_) 2 }- 1 + (A^ 2 + A w - 2 ) 



By using 



yn yn 

lim — — ^ = n, 

z->i A + — A_ 



and Eqs. fj4TT5) . fl^THjl . (14181) . fOQ"]) . we have 

limf^(z) = 1- , (1 ~ Q)(fc ' 1} , (4.20) 
z->i^ fc v ; (l-a)iV + 2a- l' V ; 

ST* 1 } ~ (l-a)N + 2a-l ' 1 j 

for any 1 < If < iV — 1. On the other hand, the well known result of symmetric classical 
random walk (that is, a = d = 1/2) and Proposition (ii) imply 

Pl N \^) = l-^. (4.22) 

Therefore, combining Eqs. fl4.20p and f!4.2ip . and Lemma [H] for a ^ 1/2 with Eq. (I4.22p for 
a = 1/2, we obtain 

THEOREM 11 For correlated quantum walk with a = d and a e (0, 1), we have 

{N) _ (l-o)(iV-fc) + (2o-l)a 
fc W ~ (l-a)N + 2a-l 

for any 1 < K < N — 1 and initial distribution ip — *[a, /?] G $. 

It should be noted that Theorem [TT1 can be also obtained by Theorem 1 in |26j . where 
a(= 1 - /3) -> a, 5i = 5 2 -> 0, pi = p 2 -> 1, 7i -> 0, 72 -> (1 - a)iV + 2a - 1, ci -> /3, c 2 -»■ a. 
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